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Abstract 

We derive the hydrodynamic equations of motion of solid and supersolid ^He, 
that describe the collective modes of these phases. In particular, the usual 
hydrodynamics is modified in such a way that it leads to the presence of a 
propagating instead of a diffusive defect mode. The former is appropriate for 
a quantum crystal and observed in recent experiments. Furthermore, we find 
that in supersolid helium there are two additional modes associated with the 
superfluid degrees of freedom. The observation of these additional modes is a 
clear experimental signature of the supersolid phase. 
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I. INTRODUCTION 



The low temperature behavior of the strongly interacting quantum liquid ^He has been 
a subject of experimental and theoretical research for decades. In 1908 Helium was first 
liquified by Kamerlingh Onnes and in 1911 he discovered a sharp maximum in the density 
at what is now commonly called the A-point.0 After that, a number of macroscopic quantum 
phenomena like superfiuid flow, second sound, the fountain effect and quantized vortices 
were observed. Phenomenological theories were developed and justified from a microscopic 
point of view.iHi Also, the famous Kosterlitz-Thouless transition was first observed in thin 
superfiuid helium films.! 

In the solid phase of ^He, which is reached only at low temperature and high pressure 
(cf. Fig. |l|), one also expects to observe macroscopic quantum phenomena because of the 
large zero point vibration of the atoms about their equilibrium position.! Because of this, 
solid helium has been termed a quantum solid. 

In such a solid, the interstitials and vacancies are effectively delocalized due to their 
ability to tunnel through the potential barriers. At low temperatures these point defects 
then form a weakly interacting Bose gas. Furthermore, the large zero point motion results in 
an unusually rapid exchange rate of nearest neighbour atoms, which may lead to large ring 
exchanges between the helium atoms.! Bose-Einsein condensation of the defects or exchange 
processes of the lattice atoms may then open two routes to a new phase of matter at low 
temperature in which long-range crystalline and superfiuid order coexist. This is called the 
supersolid phase.! 

Theoretically the existence of such a phase has since long been anticipated. However, 
it was only recently claimed to have been observed experimentally that three dimensional 
solid ^He is a spatially ordered superfiuid, or supersolid, at sufficiently low temperatures 
and densities. The experiments leading to this claim were performed by Lengua and Good- 
kind, who measured the attenuation and velocity of sound in solid ^He for relatively high 
purities and low atomic densities of the quantum crystal.0 The temperature dependence of 
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the attenuation revealed a coupling to thermally activated excitations, consistent with the 
existence of a propagating mode in the gas of point defects that is expected to be present in 
a quantum crystal.li0 Furthermore, assuming the speed of sound of the defect mode to de- 
pend on the density of defects in the same way as in a dilute Bose gas, they found a relation 
between the temperature dependence of the phase velocity and that of the defect density. 
To consistently interpret their data they then had to assume a macroscopic population of 
the zero momentum state of the point defects, i.e. a Bose-Einstein condensation of the point 
defects. Thus the phase diagram of ^He in three dimensions would be qualitatively given by 
Fig. |l|. Following a certain trajectory in this phase diagram, ^He may undergo a transition 
from the normal phase to the superfiuid phase at some temperature Tx and subsequently 
from the superfiuid to the supersolid phase at a temperature T^. As mentioned above, the 
possibility of superfiuid flow in a solid has since long been anticipated theoretically. An- 
dreev and Lifschitz were the first to attempt to derive the hydrodynamics of a supersolid 
by including the effect of Bose-Einstein condensation of the defects on the hydrodynamics 
of an ideal crystahi In addition to this pioneering work, Liu has more recently presented a 
thorough discussion of the Andreev and Lifschitz hydrodynamics.^ 

However, it was pointed out by Martin et.al. that the conventional treatment of the 
hydrodynamic equations for a classical crystal, which doesn't include defects, is neccesarily 
incomplete since it yields the wrong number of modes.0 They identified the missing mode as 
a mode in the defect density. This implies that the hydrodynamics of Andreev and Lifshits 
is also incomplete, because it does not include the non-condensed defects and as a result 
does not lead to the required defect mode in the normal state of the crystal. In addition, 
Martin et.al. assume diffusive dynamics for their defect mode. This seems to be appropriate 
for a classical but not for a quantum crystal, where the defect mode is expected to be a 
propagating mode, as is confirmed by the experiments of Lengua and Goodkind. 

Recently Stoof et.al, in respons to experiments with submonolayer superfiuid helium 
film,El derived the hydrodynamic equations for an isotropic supersolid in two dimensions 
which did include propagating behavior of the crucial defect mode.lii Moreover, the longi- 
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tudinal part of the solid hydrodynamics derived by these authors turns out to be identical 
to the system of two coupled wave equations that Lengua and Goodkind used to accurately 
model their data. However, to apply these promising results to the experiments with solid 
helium, we have to extend them in two ways. First of all we have to consider a three di- 
mensional system, and second of all we have to take into account the anisotropy of solid 
'^He, which is a hexogonally closed packed (HCP) crystal. Thus we hope to justify from 
a microscopic point of view the phenomenological equations that succesfully explained the 
propagation of sound in solid "^He and led to the first claim of a supersolid phase in this 
system. 

The paper is organized as follows. In Sec. the hydrodynamic equations describing 
a normal solid with point defects will be derived. This is achieved by deriving an action 
describing a solid with dislocations, using methods developed by Kleinert.lli From this action 
we obtain the interaction between phonons and a point defect, by seeing the point defect 
as a limiting case of a dislocation. Also, a more microscopic point of view is presented 
and dissipation is included. In Sec. fT| we then add a superfiuid degree of freedom to 
our hydrodynamic equations in the usual way and in Sec. |V] we discuss the experiment 
by Lengua and Goodkind in the light of our results. It should be noted that in order to 
understand this experiment it is not neccessary to include temperature fiucuations into our 
considerations and we will neglect them in the rest of this article. We conclude with a 
discussion and outlook in Sec. 0. 

II. HYDRODYNAMICS OF SOLIDS WITH POINT DEFECTS 

In this section we derive the hydrodynamics of a solid with point defects. This will be 
done by first considering the action describing phonons and their interaction with dislocation 
loops. We then obtain the interaction of phonons with vacancies and interstitials by shrinking 
a dislocation loop to zero radius and using a dipole-like approximation. Next, we add the 
dynamics of the point defects. The structure of the resulting theory is much like that of an 
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electron interacting with electro-magnetic fields. An intuitive microscopic picture of point 
defects is also presented which leads to an alternative derivation of the action descibing a 
solid with point defects. Finally, the hydrodynamic equations are derived. 

A. Gauge theory of phonons and dislocations 

We start with deriving a gauge theory that describes the solid phase at long wavelenghts. 
To describe point defects we only have to include dislocations into our theory. Therefore we 
can ignore higher gradient elasticity, which would be needed if we also wanted to describe 
disclinations. This section is closely related to previous work done by Kleinert0'0 but differs 
from it in the following aspects. First, we do not include higher gradient elasticity. Second, 
we consider the more general case of anisotropic solids and third, we explicitly remove the 
unphysical gauge degrees of freedom in the resulting theory of 'quantum defect dynamics'. 

The Euclidian action for a solid with dislocations of arbitrary crystaline symmetry is 
given 

S[ui] = dr J rfx - + ^{uij - ^^^' )cijki{uki - ^ (i) 

where Uij = \{diUj + djUi) is the strain tensor, Cijki is the elasticity tensor whose structure is 
determined by the specific symmetry of the crystal under consideration and p is the average 
mass density. This is the most general quadratic action compatibel with the symmetries 
of the crystal and the requirement that the Hamiltonian of the system transforms under a 
Gallilean transformation (u, t) ^ (u + vt, t) asif— s>iJ + p- v + Mv^, with p the total 
momentum of the crystal and M it's total mass. The latter determines the form of the 
kinetic energy. 

The dislocations mentioned above are topological defects, which exist because the dis- 
placement field is multivalued. In much the same way, vortices in a superfluid are a conse- 
quence of the multivaluedness of the phase field.lii The multivaluedness of a displacement 
field describing a dislocation becomes apparent when writing down what can be seen as the 
definition of a dislocation, namely 
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dui = hi 



(2) 



Here, C is a contour enclosing a dislocation line C and b is the so-called Burgers vector 
measuring the strength of the dislocation. This equation can be written in a differential 
form as follows 



where an is called the dislocation density. It is analogous to Ampere's law EijkdjBk = Ji 
and states that the displacement field is non-integrable along the line C because here the 
dislocation gives a delta function contribution. If the line C is parametrized by x(s), the 
delta function along C is defined by 



If we would use in our calculations these multivalued displacement fields, the action 
would be given by Eq. (||) with (3i = Pij = 0. However, to perform a path integral over the 
Ui it is much more convenient to use a singlevalued displacement field which takes values 
on the real axis. The unphysical singular contributions to the derivatives of a singlevalued 
displacement field which describes dislocations in a solid are compensated by substracting 
the quantities Pij and Pi. The relation between these quantities and the dislocation line C 
is conveniently visualized by the Volterra construction, which we now briefly explain. Given 
a solid without imperfections, a dislocation can be created by removing from this solid a 
volume V and drawing the boundary of the volume together, thus forming a surface S with 
boundary C, and restoring the crystaline symmetry everywhere except at this boundary C 

The singlevalued displacement field which describes a dislocation created by this con- 
struction is discontinuous across the surface S with a jump in the displacement field that is 
equal to the Burgers vector b. This discontinuity gives a delta function contribution to the 
gradient of the displacement field which is called the plastic distortion and is given by 




(3) 




(4) 




(5) 
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The integral measure is defined by dSi = SijkduX jd^Xkdudv if tlie surface S is parametrized by 
x(u, v). Furthermore, if the dislocation line C is moving with a speed v, the time derivative 
of the displacement field gives a delta function contribution of (3i = VjjSji. If we are not on a 
dislocation line, the physical values of the spatial and time derivatives of the displacement 
field should be continuous and are therefore given by 

{diu.r'^y' = d,u, - A, 

{drUjY^y' = drUj - (3j . (6) 

The value of these physical quantities equals what one would get by using the multivalued 
version of the displacement field to calculate the spatial and time derivatives-S We thus see 
that the action introduced at the beginning of this section is indeed just the classical action 
for a perfect crystal straightforwardly generalized to include dislocations.lii0 

To be able to actually calculate the interaction between the phonons and the defects we 
write the action in a canonical form. We do this by introducing two new fields by means of a 
Hubbard-Stratonovich transformation.0 Physically these field are the stress tensor aij and 
momentum density pi which are canonical to wfj'^'* = {diUjy^y^ — {djUiY^'^^ and {drUiY^^^ 
respectively. It amounts to adding to the Lagrangian density in Eq. (|l]) the quadratic terms 

2p 

and 
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Furthermore, to obtain the path integral representation of the partition function Z we 
have to add functional integrals over the momentum density pi and the stress tensor cxij. 
Note that the action containes only the symmetric part of the stress tensor, and we should 
therefore only perform the path integral over the symmetric part of a^j. Note also that c^-^; 
is symmetric under the exchanges i ^ j and k ^ I, and is defined by 
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The action we find after tliese transformations reads 

S[pi, aij, Ui] = rfr y rfx + (TijC~\iaki - ipiidrUiY'^^^ + iatju^!^^''^ . (8) 

We now integrate out the displacement field, which leads to the constraints 

drPj = diffij . (9) 

This is Newton's law. In order to automatically satisfy these constraints we rewrite the 
fields Pi and aij in terms of new fields Aij and Fij by 

= dT-Fij + EikidkAij 
Pj = diFi, . (10) 

Substituting these in the interaction, i.e. the last two terms in the right hand side of Eq. 
(1), and performing some partial integrations, we find that this part of the action can be 
written in terms of the dislocation density aij and dislocation current density Jmij = VmCtij, 
as 

S.ru[A,,F^j] =1 drj d^{-^A,a., - rF^,e^J„aj} • (11) 

In the process of rewriting the action we have ended up with to many degrees of freedom. 
These unphysical degrees of freedom manifest themselves in the fact that the new fields Aij 
and Fij are gauge fields. Indeed, the expressions for 0"^ and Pi are invariant under the gauge 
transformations 

Fij Fij + ^ikidk-^ij 

Ai, ^ Aij + {drAij - diArj) . (12) 

At first sight one might therefore think that the gauge freedom removes 12 degrees of free- 
dom. However, we note that these gauge transformations are themselves invariant under a 
gauge transformation, which reduces the number of gauge degrees of freedom. Indeed, the 
gauge transformations are invariant under 
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Kij -> Kij + diXj 

A^j A^j + drXj . (13) 



As a result the gauge freedom in Eq. (|T^) only removes 12 — 3 = 9 degrees of freedom in 
the expressions for pi and aij. Therefore we are left with 18 — 9 = 9 degrees of freedom 
in the fields Aij and Fij, which is exactly what we expect because there are 3 degrees of 
freedom present in pi, 9 in aij and the constraints in Eq. @ remove 3 of these. Note that 
we should also demand aij to be symmetric, which will remove another 3 degrees of freedom. 
This means that we end up with 6 physical degrees of freedom, corresponding to the usual 
6 phonon modes. 

In order to extract physically relevant information we will have to remove the gauge- 
degrees of freedom, i.e. fix the gauge. Before we embark on this problem however, we will 
prove the following equalities which we will need later on when deriving the hydrodynamic 
equations of motion for a solid with point defects. They are 

(uP^y) = icijiiiaki) 

{{drUiY'^y') = . (14) 

P 

The proof is given by adding to the action in Eq. source terms proportional to the 



currents Kij and Ki 



£^ dr J d^{KijuPi^''' + K,{drUiY^y'} . 

Expectation values of fiu^j"^^, {drUiY^'^^) are now easily calculated as follows 

U{ul^'\{d^niY^y^)) = /(^,^)lnZ(ir,,,ir,)|^,^.=^,=o , (15) 

where Z{Kij, Ki) denotes the partition function with non-vanishing source terms. Again 
performing a Hubbard-Stratonovich transformation we get 

S[ui] = dr J dxi^^ + ^aijcijiian - ipi {{drUiY'"^' + + i(Jij{u^j^" + c^mKh) 

~ K^.cr^MKki^ ■ (16) 



Eq. (JT^ now follow by differentiation. 

After this digression, we return to the elimination of the non-physical degrees of freedom 
present in the action S = So + Smt and reexamine Eq. (|1^). As mentioned, this gauge trans- 
formation is invariant under the transformations in Eq. (|13|). We use the latter invariance 
to choose a gauge in which A^-j = 0, which is always possible by letting Xj satisfy 



A, = -y^ A^,{T')dT'. (17) 
In this gauge our original gauge transformation reduces to 

Fij Fij + SikidkAij 

Aij ^ Aij - drAij . (18) 

In order for cTjj to be symmetric, we use part of this residual gauge freedom to choose F^j 
symmetric. In addition, we introduce the fields Xij by means of 

If we now take xin to be symmetric, cjjj will also be symmetric. In terms of these fields the 
free action becomes 



So[Fij,Xij] = dr J dx 



2p 

'^'^{(^rFij + £ikl£jmndkdmXln)Cij\i{dTFkl + Skpq^jrsdpOrXqs)^ ■ (20) 

At this point Fij and Xij both contain 6 degrees of freedom. We are thus left with 12 — 6 = 6 
non-physical degrees of freedom which somehow correspond to 6 degrees of freedom in Aij. 

To eliminate the remaining unphysical degrees of freedom, we expand the Fourier trans- 
form of the fields F^j and Xij the helicity-basis {e-J''*''}.0 If we take a direct product of 
momentum space V with itself, i.e. V the helicity basis is defined as the irreducible 
representations of the rotation group in this space. From group theory we know that they 
form a complete set .il Hence we can develop a given tensor field in this basis leading to 
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s,h 

i^..-(k) = $:eSf)(k)F(^''^)(k). (21) 



Because of the symmetry of Xij ^ij the 6 non-zero components are (s, h) — 

{(0, 0), (2, 0), (2, ±1), (2, ±2)}. To identify the surviving physical hehcity components, we 

note that the expression SikiS jmndkdmXin is symmetric, traceless and invariant under the 
following transformation 

Xln Xln + dlCn + dn^l ■ (22) 

If we choose as a basis in Fourier space {k„, e^^'^^(k), e^^'~^^(k)} and develop ^„ in terms of 
this basis, this transformation reads up to a factor k 

= ^(V2e!r^ + erf')^^ + V2er^'k^'^ + V2et'k^^^ . (23) 
From this expression we see that if we choose a new basis in which to develop xin given by 

r (2,2) (2,-2) (2,1) (2,-1) L L'| (^A) 

where 

el = ^(x/^eiJ" + eD = kk , (25) 

the components of xin corresponding to {(2, 1), (2, — 1), L'} are unphysical and disappear 
from the action because they correspond to a gauge transformation. The coordinate trans- 
formation from the old to the new basis is unitary and hence the new basis is also orthonor- 
mal and complete in the space of symmetric second rank tensors. In addition the elements 
{(2,2),(2,-2),L} satisfy 

^/eL^) = ke?/^ = hef"^ = . (26) 
11 



This means there are only 3 dynamical degrees of freedom left in diFij corresponding to the 
helicity components {(2, 1), (2, —1), L'}. 

We now Fourier transform the action and expand the fields Fij and Xtj terms of the 
basis in Eq. (^). We get 



1 

+ 2 



(2,-2) (2-2) ^L^M 



(27) 



where [ij ^ kl] denotes the part between brackets with the indicated interchange of indices. 
Up to now we have removed all but three unphysical degrees of freedom. By introducing 
and choosing a gauge in which F^j and Xij ^i'g both symmetric, the stress tensor was made 
symmetric. Hence we were left with 12 degrees of freedom. Then we identified 3 unphysical 
gauge degrees of freedom in Xij corresponding to the helicity components {(2, 1), (2, —1), L'}, 
which reduced the remaining number degrees of freedom to 9. Hence we expect a residual 
gauge freedom to be present in the above action corresponding to three unphysical degrees 
of freedom. As is apparent in the expression for 5*0, this is indeed the case and the remaining 
gauge freedom corresponds to 



(28) 



for (s, h) = {(2, 2), (2, —2), L}. We will see below that this gauge freedom is also present in 
Sint- To remove this remaining freedom we introduce three new, invariant, fields 

5 F(2,2) 



^'(2,2) ^ ^(2,2) ^ 



^'(2,-2) ^ ^(2,-2) + 



A;2 



(29) 



^2 



If we furthermore realize that 



12 



/c,e,, - ^e,. 
k 
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(30) 



we obtain the final expression for 5*0 wliicli contains precisely 6 dynamical degrees of freedom. 

Before explicitly writing down 5*0 we introduce a new compact notation which also sim- 
plifies the algebraic manipulations involved in the remainder of this article. We define the 
following quantities 



ir(2-i) 



( y(2,2) ^ 



X 



'(2,-2) 



V 



-X 



IL 



-(1) 
P ■ ■ 



I 



and 



^\xv — y^ij JlJiHjklK^kl )f 



( p(2,l) \ 

.(2.-1) 



I 



( 





^13 



\ 



p(2>2) 

,(2,-2) 



pL 



(31) 



A' 



1 
1 

2 



(32) 



Indices refering to the abstract vector space introduced above, are denoted by Greek symbols 
to distinguish them from their real space counterparts. This allows us to write Sq in the 
following way 

\ 



dk 



{2nf 









V 



k^C 



J 



(33) 



[X') 



Note the minus signs, which arise from partial integration. 

— * 

Next, we also have to write the interaction in terms of the physical fields F and The 
interaction in terms of Fij and Xij is given by 



(34) 



It is not immediately obvious from this equation that the interaction can be rewritten in 
terms of F and Therefore we will show this explicitly. After partially integrating the first 
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part of Eq. (|34D, the field Xij interacts with ejmndmOiij = £jmndm£ikidkPij which is symmetric 
and traceless and therefore has only the helicity components {(2,2), (2, — 2),L'}. Next we 
rewrite the second term term —iFijSimiJmij as follows 



_ didi\ / _ 



^Imn'Jmnj ~r I On „ I 02 ^Imn^ mnk 



(35) 



The second and third term on the right hand side of the above equation will give an inter- 
action with and F^^') as is obvious from the fact that the contractions diFi^ 
and djFij annihilate the components {(2, 2), (2, — 2), L}. The first term together with the 
interaction term involving Xij will reduce to an interaction with the new fields introduced 
in Eq. (P^D, i.e. x'- To see this we use 

(9^(5jj — didj = Simidm^ikjdk , (36) 



and substite this in the first term on the right-hand side of Eq. (pS]). We obtain 

F. 



ijd^ [e^^dpSt^.d,^^ ekstdse,rldreimnJmr 



yqk 



where we have used 



^ j (y^ qpi^P^tvjdv ^ ^kstdsd^Oiq 

i J dy. (^SqpidpEtvjdy^^^^^ ekstds^qk , (37) 

^ qrldr^lmnJmkl (^r^^qm^rn ^qn^rm)^mC^nk 

= drOgk ■ (38) 

We see that there is indeed only an interaction with the (2, 2), (2, —2) and L components of 
Fij- Together with the expansion for —iej these precisely form the fields x'^"*' ^ 

Inserting all this into Sint-, the interaction between the dislocations and the physical fields 
describing the phonon modes finally becomes 
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(eir^F(^'^) + e^'-'^F^'-^ - eSf J^,,} . (39) 

Prom this interaction between phonons and a dislocation loop we now want to derive 
the interaction between phonons and a point defect. To do so, we consider a point defect 
to be a dislocation loop that shrinks to zero radius. The Volterra construction shows that 
one creates a dislocation by removing a volume V from the crystal. When shrinking the 
dislocation loop, this volume finally ends up being the volume of a single atom. In this 
way one can thus remove or add the volume of a single atom, which results in creating a 
vacancy or an interstitial. Alternatively, one can take the long wavelength limit, in which 
the dislocation will effectively look like a point defect. However, naively applying the above 
procedure yields zero interaction between the point defects and the phonons because a point 
defect has no Burgers vector. Indeed, assuming that the wavelength of the fluctuations of 
the phonon fields are much larger than the radius of a dislocation allows for the action to 

— * 

be coarse grained. The physical fields F and x then interact with something proportional 
to 




However, the above quantity is zero and therefore there is no interaction with the phonon 
fields in a first approximation. We thus conclude that in order to find an interaction we 

— * 

need to have a gradient of the fields F and x over the radius of a dislocation loop. This is 
analogous to the well known multipole expansion in electrodynamics, where a dipole has no 
netto charge but interacts with the gradient of the electro-magnetic potentials. Assuming 
the physical fields to be slowly varying over the region of nonzero dislocation density we can 
perform a gradient expansion that leads to an interaction of the following form 

Sint drj ^ {k'a-x'Nl + kF ■ M ■ J*} , (40) 
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containing an extra factor of k. The defect density is denoted by A^a, J = {.J^^\ J^^\ J^~^^) 
are the hehcity components of the defect current density J and M^,^ is a matrix. If we have 
n defects located at {x*^")}, the defect density and defect current density Jj are given by 

n 

J.(x) = J2 ia^x(")g(")(5(x - x(")) . (41) 

n 

In general there can be both vacancies and interstitials present, and the charge q distinguishes 
between vacancies (g = —1) and interstitials (g = 1). Thus, the netto defect density A^a is 
in fact the difference between the interstitial density and the vacancy density A^i"* - A^a 
which is conserved because defects and interstitials can locally only be created in pairs. The 
defect density and defect current density therefore satisfy a continuity equation 

d^N^ = td,J, . (42) 

Up to this point, no specific crystaline symmetry has been assumed. However, the explicit 
form of the interaction is determined by the symmetry of the crystal under consideration. 
This symmetry constrains the coefficients a and the matrix M. In principle it should be 
possible to expicitly take the limit of a dislocation loop shrinking to zero in the interaction 
given by Eq. (|39|) . The symmetry of the crystal is then contained in aij and Jijk, because 
the Burgers vector can only be a lattice vector. There are however subtilities involved in 
doing this, and the vector a and matrix M will therefore be determined from symmetry con- 
siderations. Because we are especially interested in the behavior of solid ^He, we consider 
from this point on the special case of a hexagonally close packed (HCP) crystal structure. 
The associated symmetry group is Cqh, which contains rotations about the c-axis and reflec- 
tions in the afe-plane. The elasticity tensor Cijki for this symmetry containes 5 constants, the 
analogues of the Lame constants A and fi in the isotropic case. 

In what follows we need the field equations for F and %' which follow from the complete 
action Sq + Smt and are given by 

d^F = ^P-(^A'-F + i^M ■ /) - iQ^ ■ ddrN^ 
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(43) 



where 



P 

Q R 



I 



\ 



V 



(44) 



A 5t 
B C 

To determine the form of the interaction between the phonons and the defects we calculate 
the stress tensor aij resulting from a single point defect. We then note that from the 
symmetry of the crystal and the fact that a single point-defect has no orientation it follows 
that aij has to be invariant under the symmetry operations of Cgh, i.e. rotations around the 
c-axis and reflections in the afe-plane. Denoting a particular symmetry operation by 5, we 
get 



SikSjiakiiS ^x) = crjj(x) 



(45) 



which in Fourier language reads 



SikSji(Tki{S ^k) = (Jij(k) . 



(46) 



In the case of a static defect we can solve for aij, using the field equations for x'- These can 
be derived from Eq. (|^) and read 



X 



k' 



:C-^ ■ aNA . 



(47) 



With the notation introduced before Xti ^i^nd a,y are vectors, and C^^i, is a matrix given by 
Eq. (^2|) . Inserting this into the equation for a^j with -F = 0, i.e. considering a static defect, 
we get 



; 2 -i2) ^1 

aij = k el/ ■ X 



This means that Eq. (^q) translates into 

. {[5.,S,4?(5-'k)] . [C-\S-'k)] - [eif (k)] . [C-\k)]} • a = 
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(48) 



(49) 



Next we are going to translate this equation into a restriction on the coefficients a in the 
interaction. In order to do so, we must choose a particular form for the thus far unspecified 
helicity-basis. On each point of the unit sphere we choose three orthonormal vectors. One 
in the radial direction, i.e. k, the other two in such a way that the vector fields we get in 
this way are invariant under rotations about the c-axis. This cannot be done for the entire 
sphere and the points {ka,kb,kc) = {(0, 0, 1), (0, 0, — 1)} are excluded. Therefore the only 
points where Eq. (^91) is not, by construction, automatically satisfied for rotations about 
the c-axis, are indeed these two points. We only treat {ka,kb,kc) = (0,0,1), because the 
other point does not lead to any additional restrictions. In this point we choose 



.(1) 



.(2) 



.(3) 



(1,0,0) 
(0,1,0) 
(0,0,1) 



(50) 



This means that here our helicity-basis becomes 



,(2,2) 



,(2,-2) 



1 
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1 
1 




^3 



1 i 
i -1 




1 -i 
-i -1 




(51) 



Each of these three matrices transforms according to an irreducible representation of Cgh- 
Thus, according to Schur's lemma,ii they do not mix under the operator Cijki, which means 
that C^jy = if fi ^ u. We notice that these matrices transform under rotations over an 
angle a about the c-axis (we are only considering the point (0,0,1)) as 



,(2.2), 



eS(k) 



,(2,-2) 



:k) 



i2a (2,-2) 



(k) 



It now follows that 



Thus for infinitesimal rotations Eq. (|^) reduces to 



/ 



(2,2) 



i2aeij 



\ 



V 







V 



^33^ 



/ ^(2,2) \ 



,(2,-2) 



V 



. 



This equation has to be valid for all values of i and j. Therefore we get 



a(2.2)Crii-a(2.-2)c-2i = 0forz = j 



The only solution to these equations is 



^(2,2) ^ ^(2,-2) ^ Q 



which means that our interaction in first instance reduces to 



S,nt[F, x'] = 1 drj ^ {k'a^'^N*^ + kF ■ M ■ J*} . 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



Next we must determine the form of the matrix M. This is done by demanding the 
following equality to be valid 



ipdr{{diUiY^y') = ipdi{{drUi)^''y') 



(58) 



If the displacement field is singlevalued and continuous everywhere, the left and right hand 
side of this equation are two equivalent expressions for —idrSp. Therefore it should be valid 
when there are only point defects present. However, when there are dislocations present Eq. 
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(^) can be shown to be false. It is therefore not at all obvious that the equality is satisfied 
at this point, because we started out with a description including dislocations. Indeed, the 
above requirement actually gives a constraint on M, as we will see below. We first calculate 
ipdriidiUiy^y^) and find 

dk 



(2vr) 

= "/ -^/'"''{k'{F''')-p{^^'dANA) + kM,,{J,))} , (59) 

where we used the equations of motion in Eq. (|43|) , the completeness relation for the helicity 
basis in writing 



^ij V ^ ^ij ^ — y^ij ^kl ^ ^ij ^kl ) ^klmn^ 



^2) 
mn 



^1) . P + ^2) _ ^ _ /^l) _ M) , ^2) _ -<2) 

and the fact that only ef^- and ef^' are traceless. If we compare the result in Eq. (p9|) with 
the expression for ipdi{{drUi)^'^^^) , which reads 

dk 



J £vse--^'(F^) , (61) 



we see that dr{{diU,)P^y') = d^{{drU^)P''y') if 

MijJj = V2a^kJ^' 

= -V2a^drNA . (62) 

Furthermore, we take the interaction with the transverse part of the defect current density 
to be zero. This is justified by noting that the transverse part of the defect density is not 
a hydrodynamic variable. We will come back to this point in more detail below. The total 
interaction is now uniquely defined in terms of one parameter and given by 
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S^nAF. X'] = ^ / + V2drF^')N^] . (63) 

We have now completely specified the interaction of the point defects with the phonon 
field. However, only the phonon field has dynamics up to now. Because of their interaction 
with the phonons the point defects would of course effectively acquire dynamics, but one 
also expects the point defects to behave as dynamical particles if one could freeze out the 
phonon field. Roughly speaking, they would behave as particles in a periodic potential that 
could tunnel from one minimum to another. Therefore we have to add a dynamical term for 
the point defects. The most general dynamical term that describes propagating behavior of 
the defects is 

5o[{x(")}] = -\ r drY.xt\m.,dl)xf^ . (64) 

The form of the anisotropic mass m^j is constrained by the symmetry of the crystal. The 
field equations are now found by varying with respect to F, x and the positions x^") of the 
defects, and are given by 

dlF = ^p.(^A'-F- i^^ddrNj^ - iQ^ ■ adrNA 
drX = ■ {^A'F - t^V2adrNA^ - ■ dd^N^ (65) 

The equations of motion for x*^") can also be written as 

a^xf ^ = -iq''''^mi/-djBkiaM\^=^(n) , (66) 

where 

This way of writing it will prove useful when deriving the hydrodynamic equations of motion 



in Sec. Ill U 
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B. Microscopic picture 



An alternative approach to the derivation of an action which describes the coupled dy- 
namics of the point defects and the phonons is to start from a microscopic action. It describes 
the atoms constituting the crystal by their positions 

y» = n» + u(^) (68) 

relative to the sites {n*-*^} of an ideal reference lattice and assumes an isotropic, short 
range interaction V^(|y*-*^ — y*--'-'!) between the individual atoms. In this approach, it is clear 
that the hydrodynamic momentum density is g = ipdrU and, as we will see below, what 
approximations we implicitly made when we wrote down the free action of a point defect in 
Eq. (0). In first instance, the microscopic action reads 

^[u] = /^''rfr|l^p(9.u«)' + x:niy^^^-y^^'^l)| • (69) 

To explicitly include the point defects, we then decompose the displacement field into a part 
describing the phonons and a part describing the defects 

where the defects are located at the positions {x'^")(r)}. Inserting this decomposition of the 
displacement field into the action we get 




Since the positions n*^*^ correspond to the equilibrium positions of the crystal, the total 
potential ^({x^}) = Y.i<j V{\^^^ - x(^)|) satisfies 

5r({n« + u«}) = r({n«}) + C[(u«)2] . (72) 

For slowly varying displacements, the quadratic terms equals (1/2) / duMijCijkiUki- However, 
we cannot use this long wavelength result to find the interaction between the phonons and 
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the defects because defects cause fluctuations in u on the scale of a few lattice spacings. To 
proceed we therefore assume that u'^^-^ only changes significantly over a distance which is 
much smaller than the typical wavelength of a phonon. Moreover, in the continuum limit 
we can always write for the displacement of a static defect 

«f ^(y) = dJiy) , (73) 

because a defect is defined by a non-zero value of fo = / dx.Uii, i.e. the volume that is 
removed from or added to the crystal due to the presence of a defect. The function / thus 
satisfies / dyd^f = qvq, where q is either —1 or +1 depending on whether we are dealing 
with an interstitial or a vacancy. Note that in the isotropic case /(y — x) oc |y — x|~^, where 
X is again the location of the defect. Furthermore, to first order in the velocity we can write 
for a moving point defect 

drut^iy, t) ^ drut^iy - x(r)) = -drx,d,utf . (74) 



Expanding the action in Eq. (|7ll) up to second order in the displacements and making use 
of the above results the effective action describing phonons and point defects is found to be 

S[Ui] = ^ J rfx |^(9^uf )2 + ^U^ijCijklul'l^ + 

E / ^/r |-^xf)(5.m,,5.)xf) +g(")^oS.,<(x("\r) 
+g(")t;opM")]9,9,[9.< (x("), r)] + E,} , (75) 

where we have neglected contributions with m ^ n, assuming the defects to be sufficiently 
far apart to interact only through the phonon field. Furthermore, Ec denotes the energy 
associated with the creation of a defect. The microscopic action gives certain relations 
between the coefficients in this action. However, renormalisation changes these coefficients 
and we believe that it does not preserve the relations between them. Therefore they have to 
be treated as independent. This is important when trying to establish a connection with the 
action in terms of the stress tensor, as found in the previous section, which can be achieved 
by means of two Hubbard-Stratonovich transformations and following the same route as 
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before by introducing the gauge fields. Tlie result turns out to be identical and shows in 
particular that there is indeed only an interaction between the phonons and the longitudinal 
part of the defect current density. 



C. Hydrodynamics 

We can now derive the hydrodynamic equations for a crystal with point defects. The 
number of hydrodynamic modes is fundamentally related to the number of conserved quan- 
tities and the number of broken symmetries. The conserved quantities are the total mass, 
the total momentum and the netto number of defects A^a, the difference between the number 
of interstitials and vacancies. The associated conservation laws result in 5 hydrodynamic 
modes. In principle we also need to take into account energy conservation, which would yield 
an additional thermal diffusion mode.lll However, for our purposes it is relatively unimpor- 
tant and we will not consider it here. Note however that we can obtain the Hamiltonian 
from the action and therefore in principle also include this mode into our considerations. In 
addition to the conservation laws, translational symmetry is spontaneously broken, which 
results in 3 hydrodynamic Goldstone modes. Hence we expect to find a total of 8 hydro- 
dynamic modes. To find the equations of motion describing these modes we first identify 
the hydrodynamic momentum density gi with (pj) = ip{{drUiy^y^) , which is obvious from a 
microscopic point of view because locally it is just the momentum of the particles of mass 
m situated on the latice sites. It is important to note that it includes the momentum of 
the point defects, because we have constructed the physical quantity {drUiY^'^'^ that way. 
In the remainder of this article we do not explicitly include the averaging brackets, since it 
unneccessarily complicates the notation. 

We can immediatly write down the following equality 

drQi = djOij . (76) 
This equation is nothing but the constraints found in Sec. [11 A| in Eq. (^. In a perfect 
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crystal without defects, the hydrodynamic modes are the phonon modes, and their equation 
of motion is found by taking the time derivative of the above equation 

8% = d.drOi, . (77) 

Therefore we need to know drCTij which is easily calculated as 

= / J^s^'"'"' {^c.,^,etl) -A'.F- ^c,,,Kei? ■ V2a + " S)d^N^'j 



(27r)3 [ p 

= -c,. - J.) . (78) 

where we used the equations of motion for F and x'; the continuity equation drN^ = idiJi, 
the expression for Bij in Eq. (|67|) and the following equality 

^C,..e-t? .A'.F = k%,rk (^^Fi^'^^ + ^F^^'-^) + e^F'' 

k CijklkkkmFml 
^^k^^ijkl^^mFml 

= -ikkCijkigi ■ (79) 

Thus we find for gi the folowing equation of motion, describing the phonon modes and their 
interaction with the point defects 

a?* = -c«„a,(a.^-B.a„j„) . (so) 

We can spht the above equation into three continuity equations as follows 

drSp = idiQi 

dr'di = iSijkdjQk (81) 
dr9i = icijkid, |— + 1 - BkiN^ \ . 
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Note that it follows from Eq. (0), combined with the continuity equation for gi, that aij 
can be written as 

dk (dl5p £lmndm^n\ at /r,^\ 

aij = -iCijki— {-Q^ H I - ^BmNa . (82) 

To find the total of 8 modes, instead of the 6 given by the equations above, we need to 
include point defects into our hydrodynamic equations, as was first noted by Martin et.a/..0 
Thus, we have an additional hydrodynamic variable, the netto defect density A^a- Note that 
the transverse part of the defect current density is not a hydrodynamic variable because the 
momentum of the defects is not conserved. We now want to write down the equation of 
motion for A^a = Y.nQ^^'^^i'^ ~ x*^"-*). Using the equations of motion for the point defects, 
we get 

n n 

= E ^5.5/(x - x(")) + trnT^^Bkiakim^ - x^"))} . (83) 



To find the hydrodynamics of the netto defect density A^a "we have to average Eq. (|83D 
over the initial conditions. In an isotropic gas in the absence of external forces, the term 
proportional to J^n^i^^^^j^^^i^ ~ x*-"^) would be the only term present, the average of which 
is just the pressure tensor Hij. A closed set of equations giving the linearized hydrodynamics 
would then be found by writing down a gradient expansion for vTjj in terms of the hydrody- 
namic variables. In our case to lowest order the pressure tensor can only be a function of 
6Na = A^A ~ -^^A^a) because it is the only variable which is even under time reversal, and we 
get Tiij = Mij6NA + C(5A^a)- Neglecting the terms quadratic in the fiuctuations we find for 
5A^A the following linearized equation of motion 

a^A^A = -Mi.d.djN^ + im-/didjBkiaki , (84) 

where B^i B^i^N/^. Note that if we had naively introduced the dynamics of the defects 
into our theory by adding a Lagrangian density for the defect density instead of the defects 
C = —^NAilMijdidjj'^d"^ + 1)Na as was done by Stoof ei.a/.,lii we would have obtained the 
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same equation by varying the action with respect to A^a- However, our approach is more 
fundamental and illumates clearly the underlying physics of this Lagrangian density. The 
hydrodynamic equations are usually given as a set of continuity equations, i.e. with first 
order time derivatives. Therefore we rewrite Eq. (Q) as a pair of continuity equations 

drN^ = idiJi 

drJi = l^^ ( — |r^^A - ^^gl — ^^^^^^ j ■ (^^) 

We stress that this is actually only an equation for the longitudinal part of the defect current, 
= didj/d'^Jj. The transverse part is not a hydrodynamic variable and is anticipated to 
relax to zero on a microscopic time scale. This completes our discussion of the dissipationless 
hydrodynamic equations. We have obtained a set of hydrodynamic equations describing 
phonons, point defects and their interaction for a HCP crystal. They are given by Eqs. (|8TD 
and (H). 

It is interesting to note that these equations can also be derived from a hydrodynamic 
action of which the Lagrangian density is given by 

using Mijdidj = Ecm~j^didj. From the associated action, the hydrodynamic equations de- 
scribing the phonon modes and their coupling to A^a are found by writing down the field 
equations for aij and defining gi by the constraint drQi = djaij. Note that in Eq. (|8^) the 
term quadratic in aij is just /2p + aijC~jl.iaki, which is the free part of the action in Eq. (H) 
with the substition pi — > djCTij/dr- Therefore, Eq. (^61) is the analogue of the hydrodynamic 
action describing density fluctuations in a normal fluid. EE 

For completeness we write down the total set of hydrodynamic equations which as ex- 
pected amount to a total of 8 continuity equations 

drSp = idiQi 

dr^i = ieijkdjQk 
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Or9i = iCijkiOj I — I H I - BkiN/^ I (87) 

To check heuristically if we ended up with the right equations we write down the hydrody- 
namic equations in the case of a isotropic crystal and compare these to the ones found for 
a two-dimensional isotropic crystal by Stoof et.a/.El In the isotropic case Cijku Bij, and 
rriij are given by 

Cijki = >^Sij5ki + f^iSikSji + 5a5jk) 

Mij = M5ij (88) 
rriij = mSij . 

This implies the following equalities 

CijkiOjdkdi— = didp 

po^ p 

6 p 

CijkiBkidjNA = ~^^P^ \)a^diNA ■ 

Furthermore, we write the hydrodynamic equations in real time, which amount to the substi- 
tution dr —>■ —idf. As a result the hydrodynamic equations for an isotropic three-dimensional 
crystal with point defects are given by 

dtSp = -diQi 
dt'&i = —eijkdjQk 

= - ^^^di6p + ^BikidkOi - -^(/i + X)a%NA^ (90) 
dtN^ = -diJi 
dtJi = -Cpdidp - CAdiNA , 

where Cp = [ma^ /V2](4A + 4/i) and ca = M -2m{a^)^{4:\ + 3p). These are indeed the three- 
dimensional generalisations of the equations found by Stoof et.al. for the two-dimensional 
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isotropic crista! without dissipation. 

In order to give a realistic description of the system, we need to include dissipational 
effects into our hydrodynamic equations. Although there is a coupling between the phonon 
field and the defect density and thus a 'shake up' of the phonon field if a defect moves, up to 
this point there is no real dissipation because the bilineair coupling between the phonon and 
the defect modes causes mixing of these modes, but no dissipation. Therefore we include 
dissipation into our hydrodynamic equations in the standard way by first expanding the 
dissipative parti of the stress tensor to linear order in the conjugate forces and requiring the 
coefficients to be compatible with the symmetry of the system under consideration and then 
in turn expanding the conjugate forces in terms of the currents.lllBil A particularly clear 
treatment of this standard method is given by F. Jahnig and H. Schmidt.@ Quite generally, 
our hydrodynamic equations have the following form 

dt5p = -diQi 
dtOi = —eijkdjQk 

dm = -d,{a,,+a^) (91) 

dtNA = -diJi 

dtJi = -dji-Kij + vrj) , 
where the superscript D denotes the dissipative part of the 'stress' tensors, and the non- 
dissipative part has already been determined. Roughly speaking the variables 9i and Jj 
are associated with the 3 broken symmetries, whereas p, Qi, account for the conserva- 
tion of mass, momentum and defects. The most general dissipative terms allowed by the 
requirement that the time reversal symmetry of the dissipative currents is opposite to the 
associated hydrodynamic variable are given by 

= -p^ijVHdkgi + 5,,C^^'^duJk , (92) 

where we used that Ji containes only a longitudinal degree of freedom. The specific form of 
the parameters is determined by the discrete symmetries of the system, which in the case 
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of the HCP crystal ^He form the group Cgh- It should be noted that in Stoof et.al. it was 
incorrectly assumed that the transverse part of the defect current behaves as in a gas and 
diffuses to zero.El As we have seen, the correct behavior of the transverse part of the defect 
current is a relaxation to zero on a microscopic time scale. 



III. SUPERSOLID HYDRODYNAMICS 

In view of the exciting experiments by Lengua and Goodkind, our aim in writing this 
paper was also to formulate the hydrodynamic equations of supersolid ^He. Hence we have 
come to the point where we have to include into our hydrodynamic equations the superfluid 
degree of freedom. From microscopic theories developed for superfluid liquids and gases it 
is well known how we should proceed to include these additional degrees of freedom into the 
hydrodynamic equations for the normal phase.S0 First, the density p is split into a normal 
part p^- and a superfluid part p|^-, satisfying 

pS,, = pl^ + p^j . (93) 

Note that the tensorial nature of the densities is of importance in the case of an anisotropic 
HCP crystalline structure. Second, we split the total momentum density of the system into 
a normal part PijV^ and a superfluid part p^jVj according to 

= Pv? + pU^- - vj) , (94) 

where the superfluid velocity is purely longitudinal, i.e. Sijkdjvl = 0, because it is propor- 
tional to the gradient of the superfluid phase 0*. 

Furthermore, the dissipative terms have to be generalized for an anisotropic superfluid, 
and the dynamics of the superfluid velocity has to be determined. Following the standard 
treatment, the dissipative part of the stress tensor becomes 

^§ = vUi^kVi + cS^d,J, + U!ffd,pu{vt - vr) . (95) 
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The last term on the right-hand side is the most general term containing the conjugate 
variable of the phase field, i.e. dig^ = diPij{vJ — fj).!!!!!!! Furthermore, the dynamics of 
the superfluid phase field is basically determined by the Josephson relation and is given by 

dtvt = -^d,6p + P^m^ + d,C^l^d,v^ + C^'^d.djJ, + ^d,d,p%{vi - v]:) , (96) 

where Bp = p'^{dp/dp)\T,nA is the isothermal bulk modulus, p is the chemical potential per 
unit mass and (3 a = —dp/dnA\p,T- By adding the last three terms in the right-hand side of 
Eq. (p6D we have also included dissipation. However, at this point we have to realize that we 
were already dealing with a two fluid hydrodynamics in the normal solid phase, due to the 
presence of defects. This means that we also have to split the defect current density Jj into 
a normal and a superfluid part, i.e. Ji = + J^. Physically, this means that the superfluid 
current density can be caused both by the motion of defects, and by lattice vibrations.^ As 
a result we end up with the following hydrodynamic equations describing supersolid ^He 

dtSp = -digi 
dt^i = ~£ijkdjgk 

dtNA = -diJi (97) 



p 



dtvt = -^d,6p + (3Ad,NA + 
P 



+ C^'^d.d.j^ + c^'^d,d,{r^ - j;) + ^^d,d,p],{vi - vi) 

dtJ- =-^mA + (3pdi6p + 

The large number of dissipative terms makes these equations look rather intricate, but in the 
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limit — > only the non-dissipative terms remain and a considerable simplification occurs, 
as we will see below. They are easily seen to represent ten equations for the ten unknown 
quantities 6p,9i,vf,Vi, N^, Ji and Jf , realizing that ff, Jf and Jf have only one degree of 
freedom and 9i has only two degrees of freedom. 



IV. COMPARISON WITH EXPERIMENT 

We now want to compare our results with the equations used by Lengua and Goodkind 
to fit the data of their experiment in which they may have observed the supersolid phase of 
^He.0 Their phenomenological equations describe a set of two coupled harmonic oscillators. 
We show below that these equations essentially follow from our hydrodynamic equations 
describing a normal crystal with defects. 

To find the mode structure present in our dissipationless hydrodynamic equations, it is 
convenient to rewrite Eq. (|8^ in terms of the longitudinal part of the defect current density 
Jj. After taking the time derivative of the second equation of Eq. (^) and inserting the 
first equation we get 

Ji = dil — ^5 — Bkiidtaki H g^OkJk I • (98) 

To obtain a closed set of equations we then use Eq. (pS|), which expresses aij in terms of gi 
and Na- We find 



dt Ji = di 



-■DmnCmnkl —OkQl — ^klOiJi H t:;^ CfcJfc 



(99) 



We now turn to Eq. (|80D which describes the phonon modes. First we define the 
eigenvectors ^['^■'(k), n = {1, 2, 3}, of the matrix Cijkikjkk as follows 

c,,kikjkkA^'\k) = exlik)A^'\k) . (100) 

The 6 phonon modes of the ideal crystal are thus given by A\"'\'k)e'^''^^^'^^^'^^'^\ with ci;^(k) = 
A;^A^(k). In order to find the equations used by Lengua and Goodkind we first expand in 
terms of the eigenvectors A\^\'k), i.e. gi = J2n9^"'K^)Ai"'\k). We then write Ji = J^(k, T)ki 
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and insert these expressions into Eq. (^) and Eq. (^) . After contracting the first equation 
with the eigenvectors Ai and the second with ki, this leads to the following equations in 
Fourier space 

P 

dljL = -rn,^k,k, ('/3(k)J^ + l^«(")(k)(^w'j , (101) 

V P n / 

where we defined a*^"'^(k) = BijCijkikjA^"'^ and /3(k) = (Mijkikj) / {rriijkikj) — BijCijkiBki. 
Finally we consider one particular mode, say m, and elliminate the two modes with n ^ m. 
After Fourier transforming also the time variable the equations for g^'^^ with n ^ m are 
solved by 

uo^ 



Inserting this into Eq. (|101D we find 



p 



L 



-^V^ = -m,,Kk, \ /3(k) + -/-\\.)g^-^ + -^Y. f • 



P P n^ra ^ 



These equations still contain 4 separate modes. However, solutions to these equations have 
uj'^ (x k'^. Therefore we essentially find the following equations 

dlg^^^ = --k^Xlg^"^^ + fc2aM(k) 

= hm^jkj (^/5'(k) + ^!!^(^(™)^ . (104) 

These indeed describe a set of coupled harmonic oscillators and agree with the dissipationless 
limit of the equations used by Lengua and Goodkind to interprete their data. 

If we now add dissipation, the modes A*^"-* no longer diagonalize Eq. (0). However, 
there will be a new set of damped phonon modes with imaginairy eigenvalues. Proceeding as 
before, we can again elliminate two modes. We then find a coupled set of damped harmonic 
oscillators that now precisely agree with the equations used by Lengua and Goodkind. 
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To conclude this section, let us consider the dissipationless hydrodynamic equations 
describing an isotropic supersolid. The transverse phonon modes then decouple, and for the 
longitudinal part we find schematically the following equations 
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(105) 



The hydrodynamic modes can in principle be found by diagonalizing the two matrices. If we 
are in the normal phase, the first equation remains unchanged, whereas the second is absent. 
Clearly we then have four propagating sound modes. In the supersolid phase the second 
equation is also present, and we find two second sound modes in addition to the four first 
sound modes. These are however not accurately described by Eq. ( |105| ), because for that 
it is essential to include temperature fiuctuations, which we have neglected throughout this 
article. Nevertheless, it is clear from the above that to show experimentally the existence of 
a supersolid, it would be very convincing if one observes an additional resonance due to one 
of the modes associated with the superfiuid degrees of freedom. 



V. CONCLUSION 

We have derived the hydrodynamic equations for the solid and supersolid phases of 
^He. It is well known that to describe the normal solid phase, it is essential to include 
defects into the hydrodynamic equations to find the right number of modes predicted by 
the conservation laws and broken symmetries. Because we know that there are 6 phonon 
modes, the defects are usually assumed to have diffusive dynamics, giving a total of 6 + 1 = 7 
hydrodynamic modes. This is then in agreement with the 8 — 1 modes one expects from 
the usual counting argument, excluding a thermal diffusion mode. However, Lengua and 
Goodkind in their experiment observe instead propagating behavior of the defect mode. 
This brings the total number of hydrodynamic modes to 6 + 2 = 8. Therefore we introduced 
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another hydrodynamic variable, the longitudinal part of the defect momentum. We believe 
that this is justified by noting that, when counting the number of conserved quantities, 
we should also include the conservation of defects. Hence the continuity equations for A^a 
and Ji are roughly speaking associated with respectively a conservation law and a broken 
symmetry. Indeed, our equations reproduce the set of coupled wave equations which were 
used by Lengua and Goodkind to interpret their data, and lead them to the identification 
of the observed collective mode as a propagating defect mode. 

Furthermore, we have considered the hydrodynamic equations of supersolid ^He by al- 
lowing both fluctuations in the defects density and lattice vibrations to lead to superfiuid 
motion.0 If we include these superfiuid degrees of freedom into our hydrodynamic equations 
in the standard way, we end up with what one might call a four fluid hydrodynamics instead 
of the usual two fluid hydrodynamics. As a result we end up with two second sound modes 
instead of one. We expect on general grounds that including temperature fluctuations leads 
to one of these modes becoming propagating whereas the other will remain diffusive. Given 
these results it should then be possible in principle to identify experimentally an additional 
resonance in the attenuation and velocity of sound due to the coupling of these modes to 
the phonons. In our opinion this would be a more convincing experimental proof for the 
existence of a supersolid phase than the analysis made by Lengua and Goodkind. 
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FIGURES 

FIG. 1. Tentative scetch of the phase diagram of ^He. 
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